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ABSTRACT 

Surface brightness maps for the young KO dwarf AB Doradus are reconstructed from archival 
data sets for epochs spanning 1988 to 1994. By using the signal-to-noise enhancement tech- 
nique of Least-Squares Deconvolution, our results show a greatly increased resolution of spot 
features than obtained in previously published surface brightness reconstructions. These im- 
ages show that for the exception of epoch 1988.96, the starspot distributions are dominated 
by a long-lived polar cap, and short-lived low to high latitude features. The fragmented po- 
lar cap at epoch 1988.96 could indicate a change in the nature of the dynamo in the star. 
For the first time we measure differential rotation for epochs with sufficient phase coverage 
(1992.05, 1993.89, 1994.87). These measurements show variations on a timescale of at least 
one year, with the strongest surface differential rotation ever measured for AB Dor occurring 
in 1994.86. In conjunction with previous investigations, our results represent the first long- 
term analysis of the temporal evolution of differential rotation on active stars. 
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1 INTRODUCTION 

Solar-type stars exhibit signatures of magnetic activity that are as- 
sumed to be based on dynamo mechanisms operating in the star's 
outer convective zone. The detailed workings of the generation and 
amplification mechanisms of the stellar dynamos are still poorly 
understood as a result of the complex physics involved in theoret- 
ical models, and the lack of observational constraints. However, 
' it is widely accepted that differential rotation and convection are 
essential ingred ients of common theoretical amp lification models 
<ParkeJl95ilBabcoclJl9^LeightoJi963.ll969h . 

Differential rotation results from the interplay of rotation and 
convection, which leads to a redistribution of heat and angular 
momentum inside the convection zone. The thermal and density 
profiles of convective motions produce a dependence of differen- 
tial rotation on both stellar latitude and radius. It is a key pro- 
cess in the cyclic activity of the stellar dynamo being the process 
through which poloidal-toroidal field conversion occurs. To date 
there have been numerous measurements of differential rotation 
on rapidly rotating cool stars, summarised by Ba rnes et alj 12005). 
that show a steady increase in the magnitude of differential rota- 
tion towards earlier spectral types, con sistent with the theoretical 
predictions of Rtidiger *& Kiiked l2002t) . In addition, helioseismol- 
ogy reveals a differential rotation profile that varies with radius and 
latitude and is caused by the presence of an intense turbulence in 
the solar convention z one, which is driven by Reynolds stresses 
jBrun & Toomrel2002l) . 



Several methods have been used to measure differential 
rotation. The methods of IGravl <1977t) : iBrunind <198ll) and 
IReiners & SchmitJ (2002) use Fourier analysis to detect differen- 
tial rotation through line profile analysis. These methods are only 
applicable to stars that do not exhibit large cool starspots as they 
distort the shape of the line profile. Other methods reconstruct sur- 
face brightness images over a time period and then trace surface 
features to ascertain how their rotation periods are dependant on 
latitude. Exampl es include , the c ross-correlation method used by 
IDonati & Collier CameronH 19971) on AB Dor to obtain the first dif- 
ferential rotation measurement for a star other than the Sun. This 
method measures the amount of rotational shear as a function of 
latitude by cross-correlating belts of equal latitude. The sheared- 
image method extends this to include a solar-like differential rota- 
tion law into the image reconstruction process, where the rotation 
rate is allowed to vary smoothly with latitude on an image grid 
JPonati et all200d) . 

The temporal evolution of differential rotation can be deter- 
m ined by me asureme nts over severa l epochs as shown for AB Dor 
by Collier Cameron & Donati 1 2002). In this work, a matched-filter 
analysis method was used to track individual spot features in the 
trailed spectrograms. The temporal evolution of differentia l rota - 
tion on AB Dor has also been confirmed by IDonati et alj 120031) . 
through the us e of the shear ed-image method, though not for the 
same epochs as Colli er Cameron & Donati 1 2002). 

In this paper we extend the epochs for which differential ro- 
tation has been measured using the sheared-image method, by pro- 
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Epoch 



Spectral Domain (nm) Pixel Size (km s ) Resolution Telescope/Spectrograph 



1988 Dec 16,19 
1988 Dec 21 
1992 Jan 18,19,20 

1992 Dec 14 

1993 Nov 23,24,25 

1994 Nov 15,16,17 



605.9249-724.4182 
581.0809-673.5818 
498.2082-708.8381 
452.9687 -688.6151 
479.6900-751.7410 
551.0686-792.9818 



2.7 
4.6 
2.42 
2.98 
2.98 
3.56 



56000 AAT/UCL 

27000 ESO(3.6m)/CASPEC 
51100 AAT/UCL 
33600 AAT/UCL 
33600 AATAJCL 
26600 CTIO 



Table 1. Journal of observations for AB Dor showing the observed spectral domain, pixel size, resolution and the instrument used for each epoch. 



Epoch 


Julian Date 




Phase Range 




Signal-to-noise 


No of lines 


References 


1988.96 


7511.7607 


-7514.8768 


-0.37-(-0.174) : 0.381-0.672 


70-100 


80 


Collier Cameron et al. (19901 


1988.97 


7516.5617 


-7518.8339 




0.961-1.517 




280-300 


290 


Collier Cameron et al. (1990} 


1992.05 


8639.9234 


-8642.2515 


0.127-0.782 


: 0.038-0.722 


: 0.98-0.628 


50-110 


945 


Collier Cameron & Unruh (19941 


1992.95 


8970.9158 


- 8971.2737 




0.092-0.766 




110-180 


2134 


Collier Cameron (19951 


1993.89 


9314.9249 


- 9317.2631 


-0.656-0.034 


: -0.266-0.918 


: 0.175-0.886 


56-95 


1515 


Unruh et al. (19951 


1994.87 


9671.5421 


- 9673.8443 


0.087-0.591 : 


0.698-0.525 : 


0.604-0.555 


160-250 


990 


Collier Cameron et al. (19991 



Table 2. Journal of observations showing Julian dates (+2450000), phase coverage, signal-to-noise of the data set, and the number of lines used in deconvolu- 
tion. 



cessing archival AB Dor data for epochs December 1988, Jan- 
uary 1992, November 1993 and November 1994 (presented in Sec- 
tion 2). This is the first time that surface brightness images have 
been reconstructed from composite profiles computed from these 
data sets using the signal-to-noise enhancement technique Least- 
Squares Deconvolution (LSD) (presented in Sections 3 and 4). Fi- 
nally we measure differential rotation for each eopch with sufficient 
rotational phase coverage (presented in Section 5) and discuss the 
implications of our results in Section 6. 



2 OBSERVATIONS AND DATA MODELLING 

The details of the instrument configuration and observing proce- 
dures used to secure the six data sets are summarised in Table Q 
We refer the reader to the publications listed in Table|2|for further 
details. 



2.1 Data Reduction 

All frames were processed with ESpRIT, a dedicated package for 
the optimal extraction of echelle spectroscopic observations. Firstly 
a 3-D fit of the bias frame is subtracted from the flat-field and arc 
frames. Each order of the raw echelle frame is then located and 
traced using cross-correlations with a user defined reference profile. 
A linear or 2D fit to the shape of the arc lines provides the slit 
direction. Deviations from the slit direction are averaged over all 
orders provide the mean slit shape. 

The first step in the wavelength calibration procedure is to 
obtain an accurate identification of calibration lines. To start, the 
user specifies the order numbers, and approximate values for the 
wavelength of the first pixel and pixel size. This information is 
used to determine the location of lines from a calibration line list. 
A quadratic dispersion polynomial is then determined and used to 
calibrate the remaining orders. The final dispersion relation is ob- 
tained by fitting a 2D polynomial to the pixel positions and corre- 
sponding wavelengths of all lines simultaneously. The comprising 



dimensions of the polynomial fit are; 1 dimension to fit the disper- 
sion relation of each order and 1 dimension to fit its variation from 
one order to the next. 

Pixel-to-pixel sensitivity differences are removed by dividing 
each pixel in the stellar frame by the corresponding pixel in the 
flat-field. A 2D polynomial fit to the inter order background is sub- 
tracted from the stellar frame. All pixels that deviate from the av- 
erage intensity of the orde r (i.e. c osmic rays) are removed. The op- 
timal extraction of Marsh 1 1989) is then implemented. The contin- 
uum is automatically fitted firstly by using a high-degree ID poly- 
nomial, and then by a 2D polynomial to identify any systematic 
trends in the continuum's shape. A more detailed description of 
this package is given by Donati et al. (1997). 



2.2 Least-Squares Deconvolution 

LSD is a method for combining the rotation profiles of thousands 
of spectral lines in an optimally weighted manner (Donati et al. 
1997). It uses a weighted least squares algorithm to compute the 
line broadening profile which, when convolved with the known 
pattern of photospheric absorption lines in a stellar spectrum, op- 
timises the chi-squared fit to the data. The list of spect ral line s is 
obtained from the LTE model atmospheres of Kurucz 1 1993J) for 
T E //=5000 K and log g = 4.5, where features with a relative central 
depth of at least 40% of the local continuum flux are used. The total 
number of lines used for each epoch is shown in Table|2| 

It can be shown that the enhancement in signal to noise is 
equivalent to either an optimally-weighted stacking of the profiles, 
or to cross-correlation with the line pattern, scaling as the square 
root of the number of lines used. It has the advantage, however, 
that sidelobes caused by blends are automatically eliminated from 
the composite profile, giving a clean profile surrounded by flat con- 
tinuum. LSD conserves the shape of the rotational profile, implying 
that any deviations in this profile can be interpreted as brightness 
inhomogeneities on the stellar surface. 
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Figure 2. Maximum entropy fits (dashed line) to the LSD profiles (solid line) for 16+19 December 1988, observed at the AAT. The rotational phases are 
indicated to the right of each profile. 



3 RADIAL VELOCITY CORRECTION 



3.1 Telluric line alignment 

Telluric lines are used to correct for small shifts in the spectro- 
graph during the night i.e. from dewar refill and the thermal and 



mechanical relaxation of spectrograph's components, as they are 
only present in the Earth's atmosphere and are therefore at zero ra- 
dial velocity. We use the procedure of iDonati et alj J2003) in this 
analysis. Firstly the composite profile of telluric absorption lines 
in each stellar spectrum is computed using LSD, with a line mask 
comprising the wavelengths and relative strengths of known telluric 
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Figure 3. Maximum entropy fits (dashed line) to the LSD profiles (solid line) for 21 December 1988, observed at ESO. The rotational phases are indicated to 
the right of each profile. 



features (mainly of water) in the observed spectral region. The in- 
strumental velocity of the composite telluric profile is then mea- 
sured and used to define the zero-point of the velocity scale. Each 
stellar profile is then shifted by the measured instrumental velocity 
before applying heliocentric velocity corrections. 



3.2 Radial Velocity Results 

The radial velocity derived for each epoch is tabulated in Table Q 
where over the time-span of our observations there is a variation of 
3.4 km s~ 1 . The change in radial velocity reflects the orbital motion 
of AB Dor predom inantly due to the p resence of its closest com- 
panion, AB Dor C IClose et all 120051)1. AB Dor C is a lo w mass 
object (0.09 M Q ) first detected bv lGuirado J.C. et all il997l) . which 
orbits AB Dor A at a separation of 0. 156 arcseconds. In addition to 
A B Dor C, AB Dor is known to have a wi de companion AB Dor 
B Jvilhu et alll99ll:lMartin & Davevll995l) which is itself a close 
binary system, and is separated from AB Dor A by a distance of 
9.09±0.01 arcseconds. 

The AB Dor A and C orbital solution has been deter- 
mined by IClose et all l2005t) . where the parameters that define 



AB Dor A's reflex orbit are fitted by: period = 1 1.75 ± 0.25 yr, 
semi-major axis = 0.032 ±0.002", eccentricity = 0.59 ±0.03, peri- 
astron passage at 1991.8 ±0.2, inclination = 67 ±3°, fi=132±2°, 
^=107 ±7°. The orbital solution is shown in Fig.Qto fit the data 
points within the error bars and to follow the decrease in radial 
velocity at periastron (epoch 1991.8). The error is approximately 
0.5 km s _1 for the epochs of this work and 0.3 km s _1 for epochs 
after 1995, and results from broad and time variable stellar lines. 
The exceptions are the radial velocity measurements of epochs 
1992.95, 1993.89, and 1994.87, but the good fits to the line pro- 
file shown in Figures|5|&|6|for epoch 1992.05, Figure|5|for epoch 
1992.95 and Figures[lOl&[TT]for epoch 1993.89, exclude any addi- 
tional measurement errors. More empirical data points will help to 
constrain the errors of the orbital solution. 



4 SURFACE IMAGE RECONSTRUCTION 

The surface brightness images are reconst ructe d using the max- 
imum entropy code of iBrown et alJ ll99ll) and iDonati & Brownl 
( 1997). The brightness model that is incorporated in this code is 
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Figure 4. Maximum Entropy surface brightness distribution for December 1988 (epoch 1988.96), where the vertical tics at the top of the plot indicate the 
phase coverage. The plot to the left shows the fractional spot coverage per latitude bin integrated over longitude, while the plot below shows the fractional spot 
coverage per rotational phase bin, integrated over latitude. 
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Figure 1. The orbital solution of Close et al ] <2005l) for AB Dor A's reflex 
orbit. Also plotted are the empirically derived radial velocity values for AB 
Dor as determined by this analysis (epochs 1988 to 1994), The orbital solu- 
tion includes an offset in the y-axis of 3 1 .37 km s — 1 that corresponds to the 
radial velocity of the binary system, which was determined by \ 2 minimi- 
sation. The x 2 obtained for this solution is 6.023. 



the 'spot occupancy' model of Collier Cameron (1992), where each 
point on the stellar surface is quantified by the local fraction of the 
stellar surface occupied by spots. The range of spot occupancy is 
from 0, where there are no spots present, to 1, where there is maxi- 
mum spottedness. 

The imaging parameters that are used in this work for AB Dor 
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Table 3. Radial velocity measurements for each epoch of this analysis with 
previously published values for epochs after 1994. 



are stellar axial inclination, i = 60°, the projected equatorial rotation 
velocity, v sin i = 89 km s~ and the photospheric and spot tempera- 
tures respectively 5000 K and 3500 K l Don ati et all2003l). All dat a 
sets are phased according to the ephemeris of llnrrisetafl jl988); 
HJD = 2444296.575 + 0.5 1479 E. LSD profiles of slowly rotating 
standard stars (GL 176.3 and GL 367) are used as template profiles 
that represent the contribution of the photosphere and the spot to 
the shape of the intrinsic line profile. 
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Figure 5. Maximum entropy fits (dashed line) to the LSD profiles (solid line) for 18-20 January 1992, part I. The rotational phases are indicated to the right of 
each profile. 



4.1 Results 

All data sets provide good sampling of the rotational cycle of AB 
Dor. The maximum entropy images are structurally very similar 
with a polar cap and/or high latitude spots, and with varying de- 
grees of low latitude spots. The longitude resolution is approxi- 



mately 3° at the equator and the size of the smallest features in 
latitude. 

The first surface brightness image reconstructed for December 
1988 is shown in Figure|4| It shows high latitude spots that domi- 
nate over a weak polar cap and a dearth of low latitude spot cover- 
age. These sets were originally observed with the aim of detecting 
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Figure 6. Maximum entropy fits (dashed line) to the LSD profiles (solid line) for 18-20 January 1992, part II. The rotational phases are indicated to the right 
of each profile. 



circumstellar clouds on AB Dor using Ha, Ca II H and K, and Mg 
II h and k lines ICollier Cameron et alJ ll99(t) . The most striking 
feature of this data set is the fragmented polar cap. However, this 
could result from the low S/N of the data set. 

The reconstructed surface-brightness image for January 
1992 is shown in Figure The image is comparable with 



previous images of thi s data set processed without LSD 
ICollier Cameron & Unruh ( 1994) using the Ca I 643. 9nm and Ca 
67 1 .8 nm photospheric lines. It is not possible to compare exact fea- 
tures due to the signal-to-noise difference, but to a first approxima- 
tion the two images contain similar spot features at mid to low lat- 
itudes. Examples of such features are at phase 0.65, the gap in low 
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Figure 7. Maximum Entropy surface brightness distribution for January 1992 (epoch 1992.05), where the vertical tics at the top of the plot indicate the phase 
coverage. The plot to the left shows the fractional spot coverage per latitude bin integrated over longitude, while the plot below shows the fractional spot 
coverage per rotational phase bin, integrated over latitude. 
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Figure 8. Maximum Entropy surface brightness distribution for December 1992 (epoch 1992.95), where the vertical tics at the top of the plot indicate the 
phase coverage. The plot to the left shows the fractional spot coverage per latitude bin integrated over longitude, while the plot below shows the fractional spot 
coverage per rotational phase bin, integrated over latitude. 
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Figure 9. Maximum entropy fits (dashed line) to the LSD profiles (solid line) for 14 December 1992. The rotational phases are indicated to the right of each 
profile. 



to mid-latitude features at phase 0.35, and the spot features at phase 
0.5. Collier Cameron & Unruh 1 1994) also reconstructed a surface 
brightness image using the photospheric line Fe I 666. 3nm. How- 
ever, the reconstructed spot features and groupings are less simi- 
lar to the image shown in Figure which is likely to result from 
the different excitations of the two lines. At high to polar latitudes 



the images reconstructed in this analy sis show a strong and uni- 
form p olar cap whereas the images of ICollier Cameron & Unrur] 
1 1994) show weak and fragmented spot structures. It should be 
noted that the default value of spot coverage is 0.5 in the work 
of Collier Cameron & Unruh 1 1994), while in this analysis we set 
the value to be 0.999. A surface image has also been reconstructed 
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Figure 10. Maximum entropy fits (dashed line) to the LSD profiles (solid line) for 23-25 November 1993, part I. The rotational phases are indicated to the 
right of each profile. 



by Ijarvinen et alJ 120051) using photometric data for the epoch 
1991.96, which shows a primary spot at phase 0.7 and a weaker 
secondary spot at phase 1.097. While Figure^shows that there is a 
fragmented spot structure at phases 0.6 to 0.8, there is no evidence 
for a second grouping of spots at phase 1.097. 

For the single night of observations in December 1992 the re- 



constructed surface brightness i mage is broadly in agre ement with 
that previously reconstructed by Collier CameroJ J 19951) using Ca 
I 643.9 nm, Fe I 666.3 nm and Ca I 671.8 nm. Both images show 
an off-centre polar cap though the spot structure is less fragmented, 
and the lower latitude features are more clearly resolved in the sur- 
face brightness images reconstructed in this analysis. Examples of 



© 2006 RAS, MNRAS 000,[T]-?? 



Magnetic activity and differential rotation on AB Dor 1 1 



1.0 



S 0.5 



0.0 




-100 100 -100 100 
v [km/s] v [km/s] 



-100 100 -100 100 
v [km/s] v [km/s] 



Figure 11. Maximum entropy fits (dashed line) to the LSD profiles (solid line) for 23-25 November 1993, part II. The rotational phases are indicated to the 
right of each profile. 



common features include spots at phases 0.15, 0.28 and 0.45, and 
a large unspotted area at phase 0.35. For epoch 1992 .96 the photo- 
metric image reconstruction of Ijarvinen et alj 12005) shows a pri- 
mary spot at phase 0.931, and a secondary spot at phase 1.325. Due 
to missing phase coverage it is not possible to verify the position- 



ing of the primary spot from our Doppler images, but we show that 
there are no significant spot features present at phase 1.325. 

The surface brightness image for November 1993 is shown 
in Figure fl2l These images were previously reconstructed using 
the the combined lines of Ca I 643.9 nm, Fe I 666.3 nm and Ca I 
671.8 nm bv lUnruh et alj ll995l) . The images reconstructed in this 



© 2006 RAS, MNRAS 000,[T]-?? 



12 S.V. Jeffers, J.-F. Donati, A.Collier Cameron 



0.0025 0.002 0.0015 0.001 0.0005 
Frac Spot C 



CD 
CT3 



O 

O 
o 

Q_ 
CO 



0.25 
0.2 

0.15 
0.1 

0.05 




0.8 



0.6 



0.4 



0.2 



All Latitudes 





Figure 12. Maximum Entropy surface brightness distribution for November 1993 (epoch 1993.89), where the vertical tics at the top of the plot indicate the 
phase coverage. The plot to the left shows the fractional spot coverage per latitude bin integrated over longitude, while the plot below shows the fractional spot 
coverage per rotational phase bin, integrated over latitude. 



work have more intermediate and low latitude spot features. How- 
ever, there are comparable spot features at phases 0.38 to 0.55, 0.62, 
0.8 and a similar spot grouping at phase 0, though as indicated by 
the tick marks above the plot, there is no phase coverage in this re- 
gion. Similar to the comparisons of the 1992 data sets, at high and 
polar latitudes the spot stru cture is weaker and more fragmented in 
th e images re constructed by Unruh et alJ ll995T) . The surface image 
of Ijarvinen et alJ 120051) at epoch 1993.89 shows a primary spot at 
phase 1.542 and a secondary spot at phase 1.097. At phase 1.097 in 
Figure fT2l there is a large grouping of spots, while at phase 1.542 
there is a weak spot feature that is insignificant in strength com- 
pared to other reconstructed spots at phases 0.38, 0.49 and 0.62. 

The reconstructed surface-brightness image for November 
1994 is shown in Fi gurell5 | LSD has previously b een applied to the 
1994 CTIO data by Collier Cameron et alj fT999V However, in our 
image reconstructions the polar cap is stronger and less fragmented 
and at mid to low latitudes the spot features are more resolved. Ex- 
amples of similar spot features include those at phase 0.1 to 0.18, 
0. 25 to 0.35 and 0.45 to 0.55. For epoch 1994.87 the surface image 
of ljarvinen et aljl2005l) shows a primary spot at phase 1.514, and a 
secondary spot at phase 1 .097. While there is a spot feature at phase 
1.514, it is significantly weaker than its neighbouring features, and 
the spot feature at phase 0.3 in Figure fT5l We have reconstructed 
a weak spot feature at the phase 1.097, but there are other stronger 
spots close by at phase 0.3. 

At certain phases of the LSD profiles, it is possible to distin- 
guish a shallow absorption feature migrating through the line pro- 
file, which is not reproduced in the model Maximum Entropy pro- 
files. Examples of such features are shown in part II of the Novem- 
ber 1993 data set (Figure [H) at phases 0.565 to 0.815. We attribute 
these absorption features as being small regions on the stellar sur- 
face that are brighter than the surrounding photosphere in contrast 



to cooler regions that produce emission features in the line profiles. 
As the imaging code is designed only to reconstruct cool features 
on the stellar surface, it is not possible to reconstruct these bright 
features. The presence of these bright features does not interfere 
with the reliable reconstruction of cool spots. 



5 SURFACE DIFFERENTIAL ROTATION 

We used the sheared-image method of Donati et al (2000) to mea- 
sure the differential rotation of AB Dor. AB Dor is an ideal can- 
didate to measure differential rotation as its short rotation period 
(0.51479-d or 12.2053 rad d" 1 ) means that it is possible to observe 
up to two-thirds of the stellar surface in one night and to get the 
necessary overlapping phase coverage within a few days. 

The image reconstruction process also incorporates a model 
of the stellar surface whose rotation rate Q, depends on latitude ac- 
cording to the simplified solar-like differential rotation law: 

0(0) = Q eq - 5Qcos 2 9 (1) 

where Q(9) is the rotation rate at colatitude 8, f2 e9 is the equa- 
torial rotation rate and SQ is the difference between polar and equa- 
torial rotation rates. We performed a large set of of image recon- 
structions, using a two-dimensional grid of values for the param- 
eters Q, eq ,SQ. For each set of model parameters, the image recon- 
struction was driven until it reached a fixed value of spot filling 
factor. The \ 2 values of the resulting images form a "landscape" 
on this grid. The best fitting model will corresponds to the mini- 
mum in the x 2 lands cape, as models with the wrong shear will give 
poor fits to the data <Petit et all2002l) . 

The optimally fitting differential rotation parameters and their 
errors are determined by fitting the reduced \ 2 landscape grids with 
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Figure 13. Maximum entropy fits (dashed line) to the LSD profiles (solid line) for 15-17 November 1994, part I. The rotational phases are indicated to the 
right of each profile. 



a bi-dimensional paraboloid with linear and quadratic terms given 
by: 

aQ. eq + bQ, eq dQ + cdQ + dQ eq + edQ (2) 

The five coefficients are then used to solve for fl eq and dQ, 
and their errors. As discussed by Donati et al. (2003), the errors are 



determined by computing the curvature radii of the \ 2 paraboloid 
at its minimum and the correlation coefficient parameter between 
the two differential rotation parameters. 
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Figure 14. Maximum entropy fits (dashed line) to the LSD profiles (solid line) for 15-17 November 1994, part II. The rotational phases are indicated to the 
right of each profile. 



5.1 Results 

For each epoch that comprises more than one night, we converged 
the imaging code to the fixed spot filling factors obtained when 
reconstructing the surface brightness images and computed a grid 
of Q and dQ models. The resulting % 2 landscape is shown in Fig- 



ure^] for epoch 1992.05 and fitting the paraboloid (equation |2j 
to the data resulted in the values shown in Table|4| The results for 
epoch 1988.96 resulted in a map without any minimum, which is at- 
tributed to the poor quality of this data set. The temporal evolution 
of differential rotation over the epochs of this analysis is plotted 
in Figure fTTl This plot clearly shows that the magnitude of tempo- 
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Figure 15. Maximum Entropy surface brightness distribution for November 1994 (epoch 1994.87), where the vertical tics at the top of the plot indicate the 
phase coverage. The plot to the left shows the fractional spot coverage per latitude bin with integrated over longitude, while the plot below shows the fractional 
spot coverage per rotational phase bin, integrated over latitude. 



Epoch 






cos 2 9 S 




n 




(mrad d _1 ) 


(mrad d _1 ) 




(rad d" 1 ) 




1992.05 


12238.3 ±2.1 


60.1 ±5.5 


0.319 


12.219 


37761 


1993.89 


12 249.5 ±3.5 


71.1 ±9.3 


0.338 


12.226 


26645 


1994.87 


12243.1 ±3.1 


73.6 ±9.2 


0.300 


12.221 


17141 



Table 4. Summary of differential rotation parameters measured for AB Dor at each epoch, where Q. eq is the derived equatorial rotation rate at the 1<t 68% 
confid ence interval , dfl is the equatonpole differential rotation rate, column 4 is the inverse slope of the ellipsoid in the f2 eg -f2 plane (equal to cos 2 d s ref. 
IDonati etal] 120031) 1 Q s is the rotation rate at colatitude 8 S , and n is the total number of data points used in the imaging process. 



ral evolution is greater than the error bars of the differential rotation 
measurements, with the most noticeable change in dQ. of 1 1.0 mrad 
d _1 being between epochs 1992.05 and 1993.89. 



6 DISCUSSION 

6.1 Starspot distributions 

AB Dor is one of the most extensively Doppler imaged stars. In this 
analysis we have extended t he wo rk of I Donati & Collier Cameron! 
1 1997); iDonati et all Jl999h and iDonati et alj <2003l) to include 
epochs 1988.96, 1992.05, 1992.95, 1993.89 and 1994.87 to have 
a complete data set that has been processed using the same pro- 
cessing methods. One outstanding feature present at all epochs is 
a long-lived polar cap that is slightly fragmented at epoch 1988.96 
and weaker at epoch 1992.95. Direct evidence for the prese nce of a 
polar c ap has been shown by Je ffers et alj y005J) and Ueffers et alJ 
1 2006) using data from the Hubble Space Telescope of the RS CVn 



binary SV Cam. The reconstructed surface images of this analysis 
also show many small mid to low latitude spots and a spot coverage 
that ranges from 6.25% at epoch 1992.95 to 9.2% at epoch 1994.87. 
These values bear no resemblance to the results of long-term photo- 
metric measurements, which show a brightness minimum in 1988 
I Amad o"et alJ 1200 ll) . The difference between the two results can 
be accounted for by the presence of chromopheric emission, which 
for young stars such as AB Dor, show s an anti-c orrelation with the 
star's photometric variations iRadi cketalll998l) . 

For 1988.96 the reconstructed weak and fragmented polar cap 
is not conclusive due to the poor signal-to-noise of the data. How- 
ever, t here is additional evidence from the results of Kiirst er et alJ 
1 1994) that AB Dor did not posses a polar cap in 1989. This could 
show that these large polar caps could disappear periodically. If 
such behaviour is cyclic then this could be related to a change in the 
dynamo nature of the star. The weaker polar cap for epoch 1992.95 
probably results from the lower overall spot coverage at this epoch. 

The long-term stability of high, mid and low latitudes is shown 
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Figure 16. Reduced \ 2 maps showing the differential rotation profile for 
epoch 1992.05. Only points that are within 1% of the minimum x 2 are 
shown. The central black region corresponds to points where the minimum 
value of differential rotation was obtained, while the outer grey points cor- 
respond to the 1.5 cr confidence ellipse (taking each parameter separately) 
on the differential rotation parameter plain. 



to vary on a timescale that is shorter than the temporal spacing of 
our observations. T he evolution of smal l-scale magnetic features 
has been shown by iBarnes et alJ 1200 ll) to be on a time-scale of 
less than one month for the active G2 dwarf He 699. The fractional 
spot coverage per latitude bin for each epoch of this analysis is 
shown in Figure [181 The relative fractional spot coverage for low- 
latitude features (between 0° and 50°) is 24% for 1988.96, 17% for 
1992.05, 26% for 1992.95, 30% for 1993.89 and 19% for 1994.87. 
The yearly distributions as shown in Figure fT8l show that the peak 
does not become broader with time indicating that there is no ap- 
parent migration of high latitude spots towards the equator. How- 
ever this plot does show a global evolution of the spot distribution 
with possibly the polar spot becoming weaker with time and more 
spots forming at lower latitudes (compare plots for 1993.89 and 
1998.96/1992.05). In contrast to these results. ljarvinen et alj (20051 
show, from photometric data of AB Dor at similar epochs, that the 
mean spot latitude of 47-51° for 1992.05 and 45° for 1993.89 is 
significantly lower than the values of this analysis. 

The integrated brightness distributions are plotted as a func- 
tion of longitude for each surface brightness image, as shown in 
the lower panel of Figures RJQ[8]^| and !15l These distributions 
show a variation of spot coverage with longitude, and regions of 
concentrated spot coverage that could be indicative of regions of 
enhanced magnetic activity or 'active longitud es'. Active longi- 
tudes are shown to be present on AB Dor bv ljarvinen et alJ 120051) . 
where they reconstruct surface images comprising a primary and 
secondary spots from photometric data. The effect of high spot cov- 
erage can in, certain cases, result in spurious 'active longitudes' 
thoug h these are always located at the quadrature points (Jeffers 

As previously discussed we show that the Doppler images re- 
constructed here are in broad agreem ent with the location of the 
primary spot of Ijarvinen et alJ |2005), but not with the existence 
or location of a secondary spot. The difference between the results 
of this analysis and those of Ijarvinen et alJ 120051) is not a result 
of inaccurate photometric data, but due to that it is not possible 
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Figure 17. Differential rotation parameters measured for AB Dor as tabu- 
lated in Table|4|for epochs 1992.05, 1993.89 and 1994.87 (from this work) 
and 1995.94, 1996.99, 1999.00, 2000.93, 2001.99 (from Donati et al. 2003). 
The central point that is indicated by o is for epoch 1992.05 from which a 
dashed line connects the points in chronological order. The central point is 
the best-fitting differential rotation measurement. For each differential rota- 
tion measurement, the 68% confidence ellipse is also plotted. 



to correlate phases of maximum and minimum photometric bright- 
ness with regions of the highest and lowest density of low latitude 
features. This is because the shape of the lightcurve can also be 
strongly influenced by the presence of high-latitude features and an 
off-centred polar cap. The effect of including high latitude spot fea- 
tures and polar cap is shown in the plots of the fractional spot cov- 
erage as a function of longitude (lower plot of Fi gures HI 171 l8l IT2l 
and !15> where plots are shown for spot features integrated over all 
latitudes and between 0° and 50°. For epochs 1988.96 and 1993.89 
there is little difference in the general shape of the two integrated 
spot coverage distributions. However, for epoch 1992.05 there are 
notable differences at phases 0.27 to 0.4 and 0.8 to 0.95 : in par- 
ticular the degree of fractional spottedness of the plot integrated 
over all latitudes (approximately 0.02 for phases 0.27 to 0.4 and 
0.06 for phases 0.8 to 0.95), both show the same level of spotted- 
ness when integrated only between latitudes 0° to 50° when this 
is clearly not the case when integrated between 0° and 90°. The 
same traits can be seen in the fractional spot distributions of epochs 
1992. 95 and 1994.87. Sim ilar conclusions h ave also been reached 
bv lDonati et aljj2003l) and lVogt et"allh999l) . Additionally, our im- 
age reconstructions also do not show evidence for the migration of 
both primary and secondary s pots at a constant se paration or 'flip- 
flops' as shown in the work of Ijarvinen et al J l2005t> . 
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Figure 18. Variation of the fractional spot coverage per latitude bin over the 
epochs of this analysis. 



6.2 Differential rotation 

Surface differential rotation was measured for epochs with suf- 
ficient overlapping phase coverage, e.g. 1992.05, 199 3.89 and 
1994.87. In a complimentary paper, Donati et al. (2003) measure 
differential rotation on AB Dor for epochs after 1995 using the 
same method of this analysis. The resulting values of dQ and fl eq 
for all epochs are shown in Figure^] This plot clearly shows that 
the results of this paper are generally higher than previous results 
with 1994.87 having the highest differential rotation ever measured 
on AB Dor. The smaller size of the error ellipses shown in Fig- 
ure [TTlfor epochs after 1994.87 are because the data set was taken 
over a longer time period enabling a more accurate differential ro- 
tation measurement to be made. The results so far, from 1992.05 to 
2001.99, show no evidence for any cyclic behaviour. 

Another single dwarf with a similar spectral type to AB Dor 
and for which there have been multiple measurements of differen- 
tial rotation is the K2V dwarf LQ Hya. Only two measurements, a 
year apart, have been made with more than a magnitude of variation 
in Afi ranging from 0.1942 rad d" 1 to 0.01440 rad d" 1 . The large 
difference in measurements could result from the rotational period 
of LQ Hya being about three times that of AB Dor and hence the 
spatial resolution at the stellar surface is much less than that on AB 
Dor. This could result in the measurement of differential rotation 
being influenced by the short-term evolution of unresolved spots. 

The temporal evolution of differential rotation of AB Dor has 
also been determined by Colli er Cameron & Donatil 120021) . us- 
ing a method that tracks individual starspots in the dynamic spec- 
trum, for the same epochs of this analysis. The results of this 
method show a stronger differential rotation measurement for all 
epochs; 1992.05 (fi e „=12.2514 ± 0.0029 rad d" 1 and dfi=91.05 
± 13.19 mrad d" 1 ), 1993.97 (fi e „=12.2502 ± 0.0024 rad d" 1 and 
(ffi=88.49 ± 7.47 mrad d _1 ) and 1994.87 (f U=12.2481 ± 0.0047 
rad d" 1 and d«=66.84 ± 14.93 mrad d^ UDonati et alj <2003h 
discusses the discrepancy of measurements using the two methods. 



It is concluded that there are two contributing factors; (i) the value 
of vs'mi used. The spot tracking method, which is sensitive to usim, 
uses a value of 89kms~ 1 , while we use 91 kms -1 in this analy- 
sis, and (ii) the weighting of individual spots. The sheared imaged 
method of this analysis places a higher weighting on larger spots, 
while the spot tracking method places equal weighting on all spots. 
However, despite these small scale differences, the general trend of 
higher values of differential rotation for the epochs of this analysis 
are in agreement. 

Additionallv. lDonati et aljl2003t) also measure differential ro- 
tation using magnetic features, which give a different result than 
using cool spots alone. This is interpreted as being evidence that 
the dynamo is distributed throughout the convective zone and not 
confined at the base. They also compare their results with models 
of the differential rotation in the convection zone and show that the 
internal rotation velocity field is not like that of the Sun, but more 
like that of rapid rotators where the angular velocity is constant 
along cylinders aligned with the rotation axis. lDonatietalJ<2003h 
surmise that changes in differential rotation could result from un- 
derlying dynamo processes. 

The temporal evolution of differential rotation will also have 
important consequences for the stellar structure. Large variations 
will alter the spherical oblateness of the star, such that if AB Dor 
was in a close binary system it would produce long-term changes 
in the star's orbital period. This is not applicable to the AB Dor 
A/C system given their comparatively larg e (2.3 AU) separati on. 
Further support is given to the co njectu re ofjDona^taljE)03) by 
lApplegateM 19921) Lanza et alJ 1 19981) . Lanza & RodonoU 19991) 
and [Lanza 1 2005), where theoretical interpretations of the orbital 
period modulation in RS CVns is related to the operation of a hy- 
dr omagnetic dynamo in the magnetically active star. The model 
of lApplegatd fl992) assumes that these periodic modulations are 
caused by the stellar magnetic cycle converting kinetic energy in 
the conve ctive zone into larg e-sca le magnetic fields . This is ex- 
tended by Lanza et all j 19981) and Lanza & Rodonol 1 19991) to in- 
clude the effect of magnetic fields on the hydrostatic equilibrium of 
the magnetically active compon ent. These models have been fur- 
ther extended by Lanza (2005) to include an improved treatment 
of angular momentum transport in the stellar convective zone, but 
they conclude that the method of Applegate ( 1992) is not suffi- 
cient to fully explain the mechanisms of orbital migration in close 
binaries. 



7 CONCLUSIONS 

In this paper, Doppler images of the magnetically active star AB 
Dor show that its starspot coverage is dominated by a long-lived 
and stable polar cap and variable high to low latitude spot cover- 
age. The exception to this is the surface brightness reconstruction 
of epoch 1988.96 where there is evidence of a weak and fragmented 
polar cap. There is no cyclic behaviour found in either the latitude 
distribution of spots or the spot coverage fractions. Our surface 
brightness reconstructions generally show longitudes where there 
is a concentration of spots. However, we do not find a second or 
minor spot concentration which would verify the presence of ac- 
tive longitudes or flip-flop cycles on AB Dor. 

We have made the first measurements of differential rotation 
on AB Dor for epochs 1992.05, 1993.89 and 1994.87. The results 
show a temporal evolution, with epoch 1994.87 showing the high- 
est value of differential rotation ever measured on AB Dor. To a first 
order approximation the temporal evolution of differential rotation 
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show the same varia tion as the result s for th e same data recon- 
structed bv ICollier Cameron & Donati 120021) . The results of this 
work when combined with other previously published papers repre- 
sents the first long-term analysis of the detailed temporal evolution 
of differential rotation on magnetically active stars. 
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